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1 Introduction 



Hidden- variables theories, such as the pilot- wave theory of de Broglie [T] [2] 
and Bohm [3], reproduce quantum theory for a particular 'equilibrium' distri- 
bution of hidden parameters. But allowing arbitrary distributions (analogous 
to non-thermal distributions in classical physics) opens up the possibility of 
new, 'nonequilibrium' physics that lies outside the domain of quantum physics 
HISlEllTllHliniliniiniinilllllllli]. such new physics may have existed in the 
very early universe, with relaxation to quantum equilibrium having taken place 
during the violence of the big bang [U [S] [HI [Zl [HI ES] ■ In this paper, the hypoth- 
esis of early quantum nonequilibrium is shown to have observable consequences 
today. 

The concept of quantum nonequilibrium has been discussed for general (de- 
terministic) hidden- variables theories [9l [TOl [121 [SI ■ For the specific case of de 
Broglie-Bohm theory, it amounts to having configurations with a distribution P 
that differs from the usual Born-rule distribution l^"!^ (for a pure subcnscmble 
with wave function 'J) [H [SJ [51 [T3] . There were several motivations for propos- 
ing that the early universe began in a state of quantum nonequilibrium. Let us 
briefly summarise them. 

There seems to be a peculiar 'conspiracy' at the heart of modern physics, 
whereby quantum nonlocality cannot be used to send practical instantaneous 
signals. In hidden-variables theories, this conspiracy is explained as a contin- 
gency of the quantum equilibrium state. Nonlocal signalling is generally possible 
out of equilibrium (suggesting the existence of an underlying preferred foliation 
of spacetime [16]); whereas in equilibrium, nonlocal effects cancel out at the 
statistical level [51 [51 [H [TO] ■ Our inability to convert entanglement into prac- 
tical nonlocal signals is then not a law of physics, but a contingency of the 
equilibrium state. Similarly, standard uncertainty-principle limitations on mea- 
surements are also contingencies of equilibrium [S] [5] [TT] . There is a parallel 
here with the classical thermodynamic heat death: in the complete absence of 
temperature differences, it would be impossible to convert heat into work, and 
yet such a limitation would be a mere contingency of the state, and not a law 
of physics. 

Furthermore, it has been shown that relaxation towards quantum equilib- 
rium occurs, in pilot- wave dynamics, in similar fashion to thermal relaxation 
in classical dynamics (under analogous conditions and with similar caveats) 
[H [HI [SI [T71 [TB] . Given that all physical systems to which we have access have 
undergone a long and violent astrophysical history, it is then possible to under- 
stand the ubiquitous quantum noise we see around us as, in effect, a remnant 
of the big bang. 

On this view, the effectively local and indeterministic quantum physics we 
experience today emerged via relaxation processes (presumably occurring close 
to the big bang) out of a fundamentally nonlocal and deterministic physics — a 
physics whose details are currently screened off from view, by the all-pervading 
statistical noise. For as equilibrium is approached, the possibility of instan- 
taneous signalling disappears, and statistical uncertainty emerges. In effect, a 
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hidden-variables analogue of the classical heat death has actually occurred in 
our universe, explaining the above 'conspiracy'. 

The assumption of early quantum nonequilibrium was also proposed as a 
possible alternative resolution of the cosmological horizon problem (which per- 
sists even in some inflationary models [19]): the resulting early nonlocality 
might explain the otherwise puzzling homogeneity of the universe at early times 

isiaiaiiQ]. 

The search for early quantum nonequilibrium may also be motivated sim- 
ply on the grounds that de Broglie-Bohm theory (and indeed any deterministic 
hidden-variables theory) certainly allows nonequilibrium to occur. We have an 
alternative formulation of quantum physics, which yields standard quantum the- 
ory in the equilibrium limit, and which yields departures from standard quantum 
theory outside that limit. It seems natural to explore this possible new physics, 
and in particular to test for it experimentally, as far as one can. If nothing else, 
setting experimental bounds on the existence of quantum nonequilibrium can 
provide new bounds on possible deviations from quantum theory [15j . 

Finally, if hidden- variables theories are taken seriously, one is obliged to take 
the possibility of nonequilibrium seriously as well: for it is only in nonequilib- 
rium that the underlying details become visible. If the world were always and 
everywhere in quantum equilibrium, the details of de Broglie-Bohm trajectories 
(for example) would be forever shielded from experimental test. De Broglie- 
Bohm theory as a whole would then be unacceptable as a scientific theory. And 
much the same could be said for hidden- variables theories in general. 

Given the above motivations, the idea that the universe relaxed to quantum 
equilibrium from an earlier nonequilibrium state is plausible enough. However, 
to be a scientific theory it is essential to make new, quantitative predictions. 
The new physics of systems in quantum nonequilibrium has been explored in 
some detail [6j [H [9l [Til HSl [HI [S] j and a specific signature of nonequilibrium 
has been developed [HI [U . It has also been shown that if nonequilibrium were 
present at the beginning of an inflationary phase, then there would be observ- 
able consequences for the statistics of the temperature anisotropics imprinted 
on the cosmic microwave background (CMB) (TSl [10]. Further, heuristic ar- 
guments have been given, suggesting that relaxation might be suppressed for 
long-wavelength field modes on expanding space [15] (a suggestion that forms 
the starting point for the present work); and that, if relic cosmological particles 
decoupled sufficiently early, they might still be in nonequilibrium today [15] . 
However, so far, no definite quantitative predictions have been made. The aim 
of this paper is to fill this gap. 

For the first time, given a specific cosmological model, we are able to point 
to precisely where quantum nonequilibrium could be found. We accomplish 
this by studying the evolution of nonequilibrium distributions for a scalar field 
on expanding space. We show that relaxation is suppressed for field modes 
whose quantum time evolution satisfies a certain inequality. For these particu- 
lar modes, early quantum nonequilibrium is 'frozen'. For a radiation-dominated 
expansion, the inequality implies a physical wavelength larger than the (instan- 
taneous) Hubble radius. On the basis of these results, it is possible to make 
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quantitative predictions for noncquilibrium deviations from quantum theory, in 
the context of a given cosmological model. As we shall see, there are a num- 
ber of possible consequences: in particular, infra-red corrections to inflationary 
predictions for the CMB, and relic nonequilibrium particles at low energies. 

2 De Broglie-Bohm Scalar Field on Expanding 
Space 

For simplicity we consider a flat metric, 

dT^ = dt^ - a^dx^ , (1) 

where a{t) is the scale factor, H = d/a is the Hubble parameter, and is the 
Hubble radius. As is customary, we take ao = 1 today (at time to), so that |dx| 
is a comoving distance (or proper distance today). At time t, field modes have 
physical wavelengths Aphys = a{t)X, where A = 27r/fc is a comoving wavelength 
(or proper wavelength today) and k = |k| is the comoving wave number. 

Wc consider a free (minimally-coupled) massless scalar field cj) with a La- 
grangian density C = \g^/'^da(j}d°'(j) or 

C = ia302 _ \a{W(t)f . (2) 

The action is J dt J d'^x C. Wc then have a canonical momentum density tt = 
d£/d(p= a^(j) and a Hamiltonian density 

n=\''^ + \a{V<l>f . (3) 

Here, it is convenient to write the dynamics in Fourier space. Expressing ^(x) 
in terms of its Fourier components 

and writing 

for real gkr if = 1) 2), where y is a box normalisation volume, the Lagrangian 
L = J d^x jC becomes 

kr 

(For V" ^ oo, ^ X^k ^ I and VS^k- (2Tr)^6^(\i - k). Since </> is real, 
we have (/>J^ = (j>-k or gki = 9-ki, qw2 = —q-W2- A sum over physical degrees 
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of freedom should be restricted to half the values of k, for example > 0.) 
Introducing the canonical momenta 

_ dL 3. 

TTkr — -TT- = a gkr 7 

the Hamiltonian H — J d^x Ti becomes 

H = ^ i^kr , 
kr 

with ^ ^ 

Pilot-wave field theory is defined in terms of the functional Schrodinger pic- 
ture, with a preferred foliation of spacetime [SlIellllEllEalESlEllESlEe]. Here, 
the Schrodinger equation for 5" = ^'[qkr,*] is 

.9* v-^ / 1 92 1 ,2 2 \ T 



dm^ ^ d dS 



dt ^ dOkr \ dqkr 

kr 



which implies the continuity equation 



(6) 



and the de Broglie velocities 



dqkr 1 dS 



dt dqkr 



(where ^' = |5'|e*'^). The 'pilot wave' 5* is interpreted as a physical field in 
configuration space, guiding the time evolution of an individual field (/'(x, t) 
in 3-space. (Note that a similar construction may be given in any globally- 
hyperbolic spacetime, by choosing a preferred foliation [T3], so there is no need 
for spatial homogeneity.) 

Over an ensemble of field configurations guided by the same pilot wave 5*, 
there will be some (in principle arbitrary) initial distribution P[qkr,ti\, whose 
time evolution P[qkr,t] will be determined by 

'dt^^d^r\ ^d^r " ^ ^ 

kr 

If P[qkr,ti] = |^'[gkr,ti]|^, then P[gkr,t] = |*[gkr,i]|^ for all t, and one ob- 
tains empirical agreement with standard quantum field theory [31 [551 1131 [Ml HSl 
I25| . On the other hand, for an initial nonequilibrium distribution P[qkr,ii] 7^ 
|'I'['?kr, for as long as P remains in nonequilibrium the predicted statistics 
will generally differ from those of quantum field theory. In any case, whatever 
the distribution P may be (equilibrium or nonequilibrium), its time evolution 
will be given by ([7|) . 
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3 Preliminary Discussion for a Decoupled Mode 



A proper treatment of nonequilibrium freezing is given in sections 4 and 5. As 
we shall see, our treatment is applicable to arbitrary (entangled, mixed, and 
interacting) quantum states. As a preliminary exercise, in this section we shall 
discuss some elementary features for the simple case of a single decoupled mode 
k of a free field in a pure quantum state. 

From equations ([U, ([H]), and writing '5 — '!/'k(<Zki, 9k2, i)>< where >c depends 
only on degrees of freedom for modes k' k, we find that the wave function ■i/^k 
of a decoupled mode k satisfies 



V'k + T^ak (gki + '7k2) V^k , (8) 



while the de Broglie velocities for the mode amplitudes ((J'ki,'?k2) are 

^ki = , 9k2 — — (.yj 

a-^ dqki a-^ dgk2 

(with ipk = I ^k I e*'*'' ) . The time evolution of the marginal distribution pk ((Zki , 9k2 , t) 
will then be given by 

^^■^ ' ^ ' ^P.^P-)-0. (10) 



Equations ((5|)- ([TI?)) are identical to those of pilot-wave dynamics for an en- 
semble of nonrelativistic particles of time-dependent 'mass' m = moving 
in the gki — 9k2 plane in a harmonic oscillator potential with time-dependent 
angular frequency w = k/a. We may then discuss relaxation (and relaxation 
suppression) for a decoupled field mode in terms of relaxation (and relaxation 
suppression) for a nonrelativistic two-dimensional harmonic oscillator. 

Before doing so, let us recall what is already known about relaxation in 
pilot- wave dynamics. 

For a system with configuration q and wave function -0, the i7-function 

H = j dqp\n{p/\^A^) (11) 

(the relative negentropy of an arbitrary distribution p with respect to IV'I^) obeys 
a coarse-graining if -theorem similar to the classical one [H El [E] ■ Introducing a 
coarse-graining in configuration space, and assuming appropriate initial condi- 
tions for p and the coarse-grained function H[t) will begin to decrease with 
time, corresponding to an evolution of the coarse-grained density p towards 

This 'subquantum ii-theorem' formalises a simple intuitive idea: because 
p and \'ip\ obey the same continuity equation, they behave like two classical 
fluids that are 'stirred' by the same velocity field, thereby tending to become 
indistinguishable on a coarse-grained level. 
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Such relaxation has been studied numerically, on a static spacetime back- 
ground, for simple one- and two-dimensional systems [6l[8l[T7l[T8]. For an ensem- 
ble of nonrelativistic particles in a two-dimensional box, with a wave function 
consisting of a superposition of the first 16 modes, it was found that relax- 
ation occurs very efficiently, with an approximately exponential decay H{t) « 
//oe"*/*"^ of the coarse-grained i?-function (over a timescale tc) [17]. Similar 
results have been obtained for an ensemble of nonrelativistic particles in a two- 
dimensional harmonic oscillator potential [TB]. As discussed in ref. [T7], the 
numerical timescale tc was found to be in approximate agreement with a the- 
oretical relaxation timescale r defined by = ~{l/H)d?H /(M? [B]. For a 
particle of mass m, and using a sufficiently small coarse-graining length e, a 
rough order-of-magnitude estimate yields t ~ 1/ (em^/^(Ai?)'^/^) , where A£^ 
is the quantum energy spread associated with i/; |0 [17]. (The quantity t is 
analogous to the scattering time of classical kinetic theory: one expects a sig- 
nificant approach to equilibrium over timescales of order r.) If we choose a 
'natural' value e ~ 1/Ap, where Ap is the quantum momentum spread, then 
taking /^E - {Ap)^/2 m one has the simple (and rough) result 



where At is the quantum timescale over which the wave function ■0 evolves. 
3.1 Relaxation for Sub-Hubble Modes in the Minkowski 



One expects that in the short-wavelength limit, Aphys << H^^, the above 
equations ([8|)- ([T0|) will reduce to those for a decoupled mode k on Minkowski 
spacetime, because (roughly speaking) the timescale At c>c Aphys over which 
''/'k — V'k(9kij 9k2,t) evolves will be much smaller than the expansion timescale 
H-^ = a/a [11]. 

To obtain a more precise and rigorous statement, note first that at any time 
t the Hamiltonian H(t) appearing in the Schrodinger equation ([5]) has the same 
eigenfunctions and eigenvalues as are usually obtained for a two-dimensional 
harmonic oscillator of (instantaneous) mass m — a'^ and angular frequency 
oj — k/a. Thus, for quantum numbers rii, n2 = 0,1,2, ... , we have energy 
eigenfunctions 4'ni{<Iki,t)4'n2{<lk2,t) and eigenvalues -Ek(t) — (1 + ^^1 +n2)Lo{t). 
(The time dependence in 0ni(9ki,t) and 4>n2{QM2,t) comes, of course, from the 
time dependence of m = and w = k/a.) The wave function at any time t 
may then be expanded in terms of these energy eigenstates. 



r - At = l/AS , 



(12) 



Limit 



V'k(gki,gk2,t) = ^ 



C„i,n2(0'/'ni(9kl,t)(/)„2((Jk2,t) , 



"1, 12 



and the quantum energy spread Ai?k 



^{ED - {E^f will be 



AE'k = Arik • to , 
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where = ni + n2- If we consider a subsequent evolution over a time St « 
H~^, where is the timescale over which the Hamiltonian H{t) changes, 

then the Hamiltonian (together with its eigenfunctions and eigenvalues) will be 
almost constant during (t,t + 6t), and in this interval the wave function '0k will 
evolve like that of a conventional two-dimensional oscillator, with an evolution 
timescale 

_ 1 _ 1 1 _ 1 AphyS 

Ai?k Arik uj A?ik 27r 

(where we have h — 1). Significant evolution of ipk over the interval {t,t + 6t) 
can occur only if At << H^^ or 

Aphys « Ank ■ H-^ . (13) 

We may then take to be a good characterisation of the short-wavelength 
or Minkowski limit. In this limit, over timescales At = l/A_Ek << H^^, the 
wave function evolves just as it would on Minkowski spacetime. On such 
timescales, the scale factor a is approximately constant, and the equations l[8|- 
(fTO]) reduce to those of pilot-wave dynamics for an ensemble of nonrelativistic 
particles of constant mass m = a? moving in a two-dimensional harmonic os- 
cillator potential of constant angular frequency w = k/a. From the numerical 
results for the latter case [18] we may deduce that, in the Minkowski limit, for 
a decoupled mode k in a superposition |V'k) |lk) + |2k) + |3k) -I- ■•• of many 
different states of definite occupation number, the distribution Pk('i'ki, 91^2, ^) of 
the mode amplitudes will relax to equilibrium, pk ~^ IV'kl^ (on a coarse-grained 
level, again assuming appropriate initial conditions), on a timescale r given 
roughly by or 

1 1 1 

T ^ ~ . 

Aii^k A?ik LO 

3.2 Freezing of the Wave Function for Super-Hubble Modes 

In contrast, in the long- wavelength limit, 

Aphys » Ank ■ H^^ , (14) 

we have At = l/Ai^k >> H^^ and the change in the Hamiltonian H{t) over 
timescales H"^ may be treated as a sudden perturbation, leading to the con- 
clusion that the wave function '0k is approximately static — or 'frozen' — over 
timescales H^^. 

More precisely, let us again consider an evolution over an interval (t, t + 
5t) — but now with 5t of order -ff^^, so that the Hamiltonian H{t) changes 
significantly. We may write H{t -f 5t) ~ H{t) + 5H , where 6H is comparable to 
Hit). In the limit Aphys >> Ank • H^^ , the timescale At = 1/Ai?k associated 
with the 'unperturbed' Hamiltonian H{t) will be large compared to the timescale 
over which the Hamiltonian changes. We may then treat the change 5H as 
a sudden perturbation, applied over a timescale that is short compared to the 
natural timescale of the system. By standard reasoning (for example, ref. |27|). 
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we deduce that tpu hardly changes over the interval {t,t + St), that is, that ip-^ 
is essentially static over timescales H~^. 

Note that the above freezing of the wave function on timescales H^^ need 
not occur for all super-Hubble modes, since for any Aphys > H'^ the long- 
wavelength condition ([T^ will be violated if Arik is sufficiently large. On the 
other hand, of course, for any given value of Arik, the condition ([T^ will be 
satisfied for sufficiently large Aphys and the wave function will indeed be frozen. 

If ipk is frozen over timescales H~^, then the equilibrium density |V'k|^ is also 
frozen over timescales H~^. Because the evolution of |V'k(9ki, 9k2, t)\^ is driven 
by the de Broglie velocity field {qki,cik2), in accordance with the continuity 
equation 



we then expect that the trajectories {qki{t) , qk2it)) will also be frozen over 
timescales H^^. (In principle, of course, can have solutions with an es- 
sentially static density iV'kl^ and a non-negligible velocity field (9ki,'Zk2), but 
we expect these to occur only in exceptional circumstances. And in any case, 
because the phase gradient dsit^/dqur is also frozen over timescales H~^, from 
([9| we see that the velocities qkr become smaller as the scale factor a increases 
over expansion timescales .) Assuming this to be the case, it then follows 
that an arbitrary nonequilibrium distribution 7^ IV'kl^, evolving in time ac- 
cording to the same continuity equation (|15p . will also be frozen over timescales 
H^^ . In other words, at least in this simple case of a decoupled field mode, ini- 
tial quantum nonequilibrium will be frozen on timescales of order the expansion 
timescale H~^. (This is reminiscent of the well-known 'freezing' of super-Hubble 
modes in the theory of cosmological perturbations [28ll29] .) 

The above reasoning then suggests a mechanism, whereby the rapid expan- 
sion of space at early times can suppress the normal process of relaxation to 
quantum equilibrium, raising the possibility that remnants of early nonequilib- 
rium could have survived to the present day [8l I15j. However, our treatment 
so far is rather limited. We have considered only a free, decoupled mode in a 
pure quantum state. It is only expected, and not generally proven, that a frozen 
l^kl^ will be associated with a family of frozen trajectories. And, perhaps most 
seriously, while it seems significant to demonstrate nonequilibrium freezing over 
the (time-dependent) expansion timescale H^^ , in a standard — say radiation- 
dominated — expansion we have — > as i — *■ 0, so by itself nonequilibrium 
freezing over the timescale does not tell us very much about the possible 
survival of initial nonequilibrium. These limitations will be overcome in the 
following two sections. We shall first derive a rigorous condition for nonequi- 
librium freezing, applicable to an arbitrary time interval and to any (generally 
entangled) pure quantum state of a free field. Then, we shall generalise this 
condition to mixed states and to interacting fields. 



(|V'k|'<Zki)+^(|^k|'<Zk2) - 0, 



(15) 



dt a^ki 
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4 Inequality for the Freezing of Quantum Nonequi- 
librium 



To study nonequilibrium freezing over arbitrary time intervals and for arbitrary 
quantum states, we shall examine the behaviour of the trajectories themselves 
(instead of the behaviour of their guiding wave functions), thereby obtaining a 
direct constraint on the evolution of nonequilibrium distributions. 

Mathematically, as we saw in section 2, the field system is equivalent to 
a collection of non-interacting one-dimensional harmonic oscillators with posi- 
tions qicr (and with time-dependent masses m — and time-dependent angular 
frequencies cu = k/a). The Hamiltonian operator is H — X^kr -^kr, with 

TT^ 1 
TT _ 'ikr I ^„3, ,2 --2 

~ 2^ + 2" • 

Each _ffkr has (time-dependent) energy eigenvalues i?kr = ("-kr + 5)'^: where 
rikr = 0, 1, 2, (Because of the explicit time dependence in the Hamiltonian, 

the mean energy is of course not conserved: d (^h'^ /dt = l^dH /dt^ 7^ 0.) For an 

arbitrary wave functional ^'[qkr, i], the de Broglie velocity field is given by ([6]), 
and the evolution of an arbitrary ensemble distribution P[q\^rit] will be driven 
by this velocity field via the continuity equation ((7|). 

Note that the use of a classical spacetime background must break down in 
the limit t ^ Q. The equations defining our model can be trusted only down to 
some minimum initial time ti. For example, very optimistically, one might take 
the 'initial time' to be of order the Planck time, ti ^ tp ^ 10^^'^ s. 

Now, an initial nonequilibrium distribution P[<?krj^i] 7^ |^[9krj ^i] P can in 
general relax to equilibrium (on a coarse-grained level) only if the trajectories 
wander sufficiently far over the region of configuration space where is con- 
centrated; otherwise, for example, if P were initially small in regions where j^fp 
is large, P could remain so, and equilibrium would never be reached. We may 
then write a simple condition for initial nonequilibrium to be 'frozen', by consid- 
ering the displacements of the trajectories, and requiring that the (equilibrium) 
mean magnitude of the displacements be smaller than the width of the wave 
packet. 

Let us write the total configuration of the system as q{t). Note that 
is in general an entangled function of all the Qkr's. Even so, given the initial dis- 
tributions P[q, ti] and one may calculate the corresponding marginals 
for just one Qkr (for some given kr). If the resulting two marginals are equal or 
unequal, we may say that we have equilibrium or nonequilibrium respectively, 
for the given degree of freedom gkr- In this sense, it is clearly possible for some 
of the gkr's to be in nonequilibrium while the others are in equilibrium. 

Let us now consider the motion qwrit) of one degree of freedom, for some 
given kr, over a time interval [ti,tf]. An initial point qi^riU) undergoes a final 
displacement (5gkr(i/) = //.^ dt q-^rit), where the velocity gkr is given by (l6|)ll 

^Note that trajectories in one-dimensional g^r-space do move past each other, being com- 
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Let Agkr(i) be the width — with respect to gkr — of the quantum distribu- 
tion t]p at time t. If the whole family of trajectories qwr{t) (with fixed 
kr and arbitrary initial total configurations q{ti)) were such that the magni- 
tude \5qy:r(tf)\ of the final displacement were small compared to Agkr(i/), then 
relaxation (with respect to gkr) during the interval would in general be 

impossible, as the configurations would not move far enough for the two 'fluids' 
P and l^'p to be significantly 'stirred' or mixed (with respect to (/kr)- This is 
clear because the time evolutions of P and j^/p are determined by the same 
continuity equation and the same family of trajectories. For example, if j^fp 
is initially spread over an interval [a,b] of gkr-space of length ^ A(7kr(ii), and 
if the displacements of all the trajectories during [ti,tf] are indeed such that 
\Sqkr{tf)\ « A(7kr(i/), then j^ip will essentially remain spread over [a,b] dur- 
ing [ti,tf] (with A5kr(i/) ~ AgkT-(ii)); while if P is, say, initially confined to 
the left half of the interval [a, b], it will essentially remain so during [ti,tf], and 
there will be no significant evolution towards equilibrium (for the coordinate 

Thus we might take our condition to be |(5<Zkr(i/)| << ^Ikritf). How- 
ever, if there were some isolated trajectories for which \Sqkr{tf)\ ^ A(7kr(i/), or 
even \Sqkr{tf)\ ^ ^Qk.r{tf), while most trajectories still satisfied |(5gkr(i/)| << 
A(7kr(i/) (where 'most' could be defined for example with respect to the Lebesgue 
measure or with respect to the |^p-measure), then relaxation would still be im- 
possible in general. Hence we may take the weaker condition 

(l'59k.(t/)|),q« Agkr(t/) , (16) 

where {\6qkritf)\) is the average of \Sqk.ritf)\ over an equilibrium ensemble. 

The condition (fT6|) implies that 'most' of the ensemble cannot move by 
'much' more than a small fraction of A(7kr(i/), in the following precise sense. 
Define S = \6qkr{tf)\ = \qkr{tf) — qkr{ti)\ > (where S = 6{qi,tf) is a func- 
tion of the initial total configuration qi = q{ti)). From (|16p . wc can write 
< eA(7kr(i/) for some e « 1. We can then show that 'most' values 
of 5 cannot be 'much' bigger than eAgkr(i/) — where we define 'most' with 
respect to the equilibrium measure iijpdg.i over the ensemble of initial 

configurations g^, and where we define 5 to be 'much' bigger than eA(7kr('^/) if 
5 > 2eA(7kr(i/). Let R be the set of initial points qi such that S > eAqkr{tf) + d, 
for some fixed d > 0. Such points make up a certain fraction F of the ensemble, 
that is F = Jj^ dQi l^i^ii^dP- We have a mean 

(^)cq= / dq^ mq^,U]\'■5{q^,tf) . 



ponents of higher-dimensional trajectories q{t) (unlike in a strictly one-dimensional system, 
where the single- valuedness of the velocity field prevents trajectories from crossing). 
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Since 5 > for all qi, we have 

('5)eq> / dq., \^[q,,U]\\6{q,,tf) 
JR 

> / dq,\'f[q„t,]\'^.{eAqkr{tf) + d) = F.{eAqi,r{tf)+d) . 

JR 

Given (fTB|) . or (S)^^ < eAgkr (i/), we then have eAgkr(^/) > F. (eA(7kr(^/) + d) 
and so 

(1 - F) 

d<^-^r^eAqUtf) ■ (17) 

Now, if > i (that is, if 'most' initial points yield S > eAqi^r{tf) + d), then 
{1 — F) / F < 1 and so d < eAq]^r{tf)- We may then indeed conclude that 'most' 
of the initial ensemble cannot move by 'much' more than eA(7kr(i/)- In this 
case, even an approximate relaxation cannot (in general) occur. 

If (fTB]) is satisfied, then, relaxation will in general be suppressed. Of course, 
while (fT6|) is a sufficient condition for relaxation suppression, it is not neces- 
sary: in principle, the trajectories could even wander over distances larger than 
A(7kr(i/) but without a sufficiently complex flow to drive the ensemble towards 
equilibrium. (As discussed in section 7, it is reasonable to assume that this is 
unlikely.) 

While (|16p provides a condition for the freezing of quantum nonequilibrium, 
in practice it is likely to be more stringent than is necessary. Without attempting 
to give a rigorous justification, we expect that there will be cases where the 
weaker condition 

{\Sqkr{tf)\)^^<Aq^r{tf) (18) 

suffices to prevent relaxation, at least partially (that is, some significant relax- 
ation towards equilibrium will occur but significant deviations from equilibrium 
will remain). Generally speaking, we expect that the transition from essen- 
tially complete relaxation suppression to essentially full relaxation will take 
place when the ratio r = (|^gkr(i/)|)(,q /^'?kr(i/) increases from r « 1 to 
r > 1, with the critical demarcation line being somewhere in the neigbour- 
hood of r ^ 1. We therefore expect that the weaker condition (fTS]) will define 
(approximately) essentially the whole of the suppression regime, including those 
cases where significant relaxation towards equilibrium does occur but where sig- 
nificant deviations from equilibrium still remain. (Note that (llSp implies that 
'most' of the ensemble cannot move by 'much' more than A(7kr(i/), in the sense 
given above.) 

Pending a more precise treatment, then, here we shall take ^TE\\ as our condi- 
tion for the freezing — or at least partial freezing — of quantum nonequilibrium. 

Let us now proceed to draw inferences from (fT8|l . Note first that the final dis- 
placement (5(7kr(i/) has modulus |(5gkr(i/)| < It^ dt |qkr(OI (where J^l^ dt \qkr{t)\ 
is the path length). The equilibrium mean {\Sqkr{tf)\)^^ then satisfies 

{\Squr{tf)\),^ < l^f^' dt |gkr(t)|^ = f^' dt (|gkr(i)l>cq , (19) 
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where the equilibrium mean speed (|(7kr(OI)cq time t is 

(I'Zk.WDea = / dq \nq,t]\^\qUq,t)\ (20) 



'cq 

(the velocity q]^r{<l,i) being given by ^ as a time-dependent function of the 
total configuration q). 

For the sake of clarity, let us explicitly demonstrate the last equality in (fT9|) . 
The initial equilibrium distribution represents an ensemble of initial 

(total) configurations qi. From each q^, the de Broghe velocity field generates 
a trajectory q{t) (for the whole system), and each such trajectory implies a 
subsystem trajectory gkr(0- Thus, at any time t, the subsystem velocity gkr 
may be regarded as a function of qi and of t (assuming the wave functional 
is given). We may then write qkr = 'Zkr(9i,i) — where of course qkriqi,t) and 
9kr('?, t) here denote two different functions of the first argument. (This notation 
is strictly speaking ambiguous, but clear from the context.) We then have 

(l'?kr(i)|)eq = J dq^ |*fe,i,]|2|9kr(gz,t)| (21) 

(with the mean taken over the distribution at the fixed initial time 

ti), so that 

' dt (|gkr(i)|),q = j dq, {^l' dt \ql,riq^,t)\ 

dt \qkrit)\ 



as used above. (We have shifted notation back and forth, with qhr(t) and 
%Lr{(li,t) denoting the same thing.) 

Using {x) < \J (x^) for any a;, we then have 



Now note that, at any time t, 



dq] 



kr 



(where (^Q^ denotes the usual quantum expectation value for an operator Cl). 
The last equality follows from 

7kr / J ^9kr dqkr 
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and from 

Thus, since (9|^'|/i9gkr)^ > 0, we have 

a'{\q^r\'),^<{nl) , (23) 

and so 

{\Sq^r{tf)\l^< r dt^^^^) (24) 



(where it is understood that quantities under the integral sign are evaluated at 
time t). 

Since (g^r) > 0, we also have 



(tt^^) < 2a3 {H^^^ , (25) 

and so 



Introducing the number operator nkr, where (nkr) > 0, the mean energy in 
the mode kr is 

(Hkr) - ((rikr) + 



2' a 
We then have 



(l'5'?kr(t/)|),q < / dt — V2fc((nk.) + l/2) . (26) 

The mean (|<^'Zkr(^/)|)oq time is to be compared with the width Aqkr(i/) 
(with respect to gkr) of the quantum distribution |5'[(7,i/]p at time tj. Using 
the uncertainty relation AgkrAi'kr > 5 and ATTk,- < •\/ (vr^^) , we have 1/Agkr < 
^V\^kr)- Again using ([25]) we then have 



1/Aqkr < 2 J2a3 (ilk^) = 2ay/2k{{h^r) + 1/2) . (27) 



Combining the results (|26p and (|27p . we obtain an upper bound for the ratio 
{\Sq^ritf)\), 



Agkr(i/) 



< 4kafJ{h^r)f + 1/2 / ' -4\/<^kr) + 1/2 (28) 

' Jti ^ 



(where a/ = a(tf), and so on). Note that (nkr) is in general a function of time 
t, and that the inequality ((28|) holds for any arbitrary (in general entangled) 
state 
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We may now consider the following inequality, that the right-hand side of 
28|) is less than one, that is 



- > 4a/y'(rikr)j: + 1/2 j dt ^^("kr) + 1/2 



(29) 



When this 'freezing inequality' is satisfied, (|(5gkr(^/)|)oq /^'J'kr(^/) < 1 and 
initial quantum nonequilibrium will be (at least partially) 'frozen'. 
We may also write ([2^ directly in terms of ( iJ^r / , yielding 



4afJaf{Hi,r)^J^ 'dt ^\la{H^r) < 1 • (30) 



The dependence on the wave number k is of course still present in yi/kr 
Roughly speaking, the freezing inequality ([SU]) requires that the mean energy 
Hkr ) hi the mode kr be not too large over the time interval [U, tf] (see below) 



5 Generalisations 

Before discussing the consequences of the above results, let us first generalise 
them to more realistic situations. The above derivation of the freezing inequality 
(|^^ (or ||5D|) ) assumed that the quantum state was pure and that the field was 
free. The derivation is easily generalised to mixed states and to (finite models 
of) interacting fields. 

With these generalisations in hand, one can then discuss nonequilibrium 
freezing for a mixed (for example thermal) ensemble of interacting particles, 
and one can apply the results to realistic models of the early universe. 

5.1 Mixed States 

In quantum theory, a mixed state is represented by a density operator p, which 
may be written as a decomposition 

/5 = ^p„|*„)(vl/„| , (31) 

a 

with appropriate probability weights Pa and pure states For a scalar field 

4>, the quantum-theoretical distribution for cf) will be 

PQT[0,t] = (0|p(t)|<^) -^Pa|vl'„[0,i]|2 . (32) 

a 

The decomposition of p is generally non-unique, and different decompositions 
of the same p are physically equivalent in all respects. 

The situation is different in pilot-wave theory. A mixed quantum state is 
interpreted as a statistical mixture of physically-real pilot waves ^a, with prob- 
ability weights Pa, corresponding to a preferred decomposition of p [30j . For a 
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given element of the ensemble, the de Broglian velocity of the actual configura- 
tion is determined by the actual pilot wave ^q. A different decomposition of 
p would generally yield different velocities, and so be physically distinct at the 
fundamental level. (Note that, in quantum nonequilibrium, the velocities and 
trajectories for single systems can be measured without necessarily disturbing 
the wave functions dUmj, enabling the preferred decomposition to be detected. 
The operational equivalence of different decompositions of p is a peculiarity of 
the quantum equilibrium state; see ref. [13].) 

Now, given such a preferred decomposition, for each pure subensemble with 
wave functional t] — taking the system to consist of a scalar field 4> — 

we may define a distribution Pq[0, t] (generally ^ |^a[</),t]p) and an associated 
7?-function 

(for some appropriate measure Dcj)). The whole ensemble has a distribution 

P[(t>,t]=Y.PMM , (33) 

a. 

and the mean i/-function 

obeys a coarse-graining 7J-theorem (for a closed system with constant pa) [13j . 
The equilibrium minimum H — (which may be approached in a coarse-grained 
sense) corresponds to = and Pa = \'^a\'^ for every a, so that ([55)1 reduces 
to dSll). 

Thus, we may discuss relaxation for a mixed state in terms of relaxation for 
its component pure subensembles. We may then consider the freezing inequality 
([^^1) (or (|5D|) ) for each pure subensemble separately. Clearly, the inequality 
might hold for some subensembles and not for others (or for all of them, or 
none). 

If the (quantum) mean occupation number for state |5'q) is {fikr) a = {^a\ ^kr 
then for a mixed state pip the overall mean occupation number will be 

a 

For example, for a thermal ensemble with temperature T, we will have the 
Planck distribution ^ 

(n.kr)p = ^hu^/k^T „ I ■ 

In general, {fikr)^ ^oi a pure subensemble will differ from (nkr), and the total 
ensemble will contain a range of different values for (nkr)^. Initial quantum 
nonequilibrium will be frozen for the pure subensemble with wave functional 
\E'a, if the corresponding quantity {ny^r)^ satisfies the freezing inequality (|29p 
(with (jT-kr) replaced by {•hkr)^)- 
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To investigate which (if any) pure subensembles will satisfy the freezing 
inequality ([29]) . we need to know the quantities {fi^nr) a functions of time, that 
is, we need to know which pure states \^a) are present in the total ensemble. 
Despite the operational equivalence of different decompositions of p in quantum 
theory, it has been argued that, in the case of thermal (canonical) ensembles, 
there is a natural probability measure on the space of normalised wave functions, 
the 'Gaussian adjusted projected measure', which is unique for each p, and which 
may be used to define a preferred decomposition [M]. This proposal has been 
applied to the case of an ideal gas (though described in terms of particle theory 
rather than field theory) [32]. For our purposes, we would need to apply the 
preferred measure to a thermal ensemble of wave functionals in field theory on 
expanding space, and use the results to deduce which (if any) subensembles of 
finite measure satisfy the freezing inequality ((29|) (or ((30)) ). We do not attempt 
such a calculation here, but it should be clear that the problem is well-defined. 

If certain pure subensembles — with labels a in some set S — are pre- 
dicted to be frozen, then (assuming initial nonequilibrium) the total ensemble 
distribution of <j) will take the form 

where Pa ^ l^'aP (for a E S), and P[<j>,t] will generally differ from the equilib- 
rium result ([32]) . 

The physics of nonequilibrium mixed states needs further development. In 
particular, one should explore how measurements could probe the nonequilib- 
rium physics particular to a specific pure subensemble (noting again that, unlike 
in quantum theory, in nonequilibrium pilot- wave theory it is operationally mean- 
ingful to speak of the physics of component pure subensembles). However, the 
above suffices for the purposes of this paper. 

5.2 Interacting Fields 

Our derivation in section 4 of the freezing inequality ([29)1 assumed that the field 
(j) was free. The derivation is easily generalised to interacting fields, at least if 
one considers finite models with an appropriate high-frequency cutoff (so that 
divergences may be ignored). 

Let the scalar field (p interact with other fields, denoted collectively by <&. 
(These other fields need not be scalars.) We have a total Hamiltonian 

Htotai = H + H<j, + Hi , 

where H and are respectively the free Hamiltonians for (f> and <&, while Hi 
is the interaction Hamiltonian. 

We may still of course write (j) in terms of its Fourier components (jj^r, and 
the free Hamiltonian H still decomposes into a sum H = X^kr ^kr, with Tfkr = 
(rikr + 5)|i exactly as before. Equation ((22)) still holds (for a pure subensemble 
with wave functional ^P, and where the total configuration q now includes $ as 
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well as (p). So we still have the inequality ([^5]) . The other inequalities — such 
as (1251) and Aq^r^T^kr > ^ — are also valid as in the case of a free field. We 
therefore arrive again at the upper bound (|28p and the freezing inequalities (|29p 
and pHl. 

The only difference from the free case is in the time evolution of (fikr) (or 
of (^Hkr^), which now involves contributions from Hj: 

The calculation of ( function of time t will then be more complicated 

than in the free case, where only the first term appears on the right hand side. 
(The evolution of {n-k.r)t in the free case is studied in section 8.) 



6 General Implications of the Freezing Inequal- 
ity 

Quite generally, then, even for an interacting field in a mixed state, we may 
conclude that relaxation will be suppressed — that is, nonequilibrium will be 
frozen — for modes whose (time-dependent) mean occupation number (fikr) 
satisfies the inequality ([29|) . 

For a given time evolution, defined by a{t) and (nkr)^ (for all kr) on [U^tf], 
it is of course possible that (|29p will not be satisfied for any value of fc, and 
that all modes relax (at least approximately) towards equilibrium during the 
interval [ti, tf]. On the other hand if, for a given time evolution, ([29]) is satisfied 
only for certain values of fc, then we can predict that significant deviations from 
quantum equilibrium are to be expected only for those particular values of k. 

We emphasise that, for each mode, whether or not the inequality is 
satisfied depends on the history of the expansion and on the time evolution of 
the quantum state of the field. 

For a radiation-dominated expansion on [ti,tf], with a{t) = af(t/tfY^'^ , we 
may make a general statement about the kind of modes that can satisfy ([29|) : 
the physical wavelength Aphys(i/) = a/(27r/fc) at time tf must be larger than 
the Hubble radius HJ^ at time tj (assuming that tf > (1.17)i,;). 

This is easily shown for any quantum state. Since (rikr) > 0, the inequality 
((29|) (assuming it to hold) implies that 

1 /■*•' 1 2tf 

T>2af dt -^^Htf/t,) , 

k Jt^ af 

or 

Aphys(t/) > 2TrHJ^ Htf/U) , (35) 

where HJ^ = 2tf and where the right-hand side is indeed larger than Hy^ if 
tf > iiexp(l/27r) ~ (1.17)^^. (This is of course not to suggest that the freezing 
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inequality is satisfied for all super-Hubble modes: rather, if the inequality is 
satisfied, then the corresponding modes must be super-Hubble.) 

Note that, in any reasonable application of this result, the factor \n{tf/ti) 
will not be large. For example, taking ~ ~ 10"*"^ s, for ~ 10"'^^ s 
(the time at which inflation begins in some models [2S]) we have \n.{tf/ti) ~ 
In 10* ^ 20, while even for tf ^ 1 s (the time of neutrino decoupling) we have 
ln{tf/ti) - InlO''^ - 10^ The factor 2TT\ii{tf /U) is then likely to be at most of 
order 10^ — 10'^, in which case the minimal value of Aphys(i/) for nonequilibrium 
field modes will be at most two or three orders of magnitude larger than the 
Hubble radius HJ^. 

On the other hand, again for a radiation-dominated expansion, the true 
lower bound on Aphys(^/) (set by (|29|) ) will be much larger than 2nHJ^ \n{tf /ti) 
if {nkr)t >> 1 during the period [ti,tf], as is clear from ([29)) . 

Thus, de Broglie-Bohm theory (with the assumption of early quantum nonequi- 
librium at some initial time ti) predicts that residual or 'frozen' nonequilibrium 
will exist at later times tf > U for modes satisfying the inequality ([^^ , where for 
a radiation-dominated expansion the physical wavelength Aphys(i/) of nonequi- 
librium modes at time tj must be bigger than 2itHJ^ \n{tf/ti) (at least). 

If we take the freezing inequality in the form pop . we see that, roughly 

speaking, it entails an upper bound on the mean energy ( _ffkr ) per mode over 



time. More precisely if, for example, (/fkr) > (-ffkr) throughout [ti,tf 
then ([50)1 implies that 



mm 



^kr) < . (36) 

For a radiation-dominated expansion, and assuming tf/ti » 1, we then have 
(inserting h) 

'Hur) < ^ ^ ^ (37) 

\ /min IGtf 

(where, dimensionally speaking, HJ^ — 2tf is the Hubble time and cHj^ is the 
Hubble radius). 

Finally, we note that violation of the freezing inequality (|29p in the infra-red 
limit k ^ requires that (rikr) be divergent as fc ^ 0. Alternatively, for (|3Gp 

to be violated as fc — > 0, the mean energy per mode ^-ffkr^ must remain finite 

as fc — > 0. 



7 Relaxation for Modes Violating the Freezing 
Inequality 

We have shown that, for modes satisfying ((29l) . relaxation will be suppressed 
over the time interval [ti,tf]. For a radiation-dominated expansion we know 
from (|35p that such modes, if they exist, must have super-Hubble wavelengths. 
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Further, as discussed in section 3, we know from previous studies that relaxation 
is Ukely to occur in the short- wavelength (Minkowski) limit. What can we say 
about modes that violate the freezing inequality ([M)) without approaching the 
Minkowski limit? 

Our derivation of the upper bound (j28p made use of several general inequal- 
ities (such as (tt^^) < 2a^ ^-^kr^)- For a large class of quantum states, these 
general inequalities could be replaced by approximate equalities, to be used as 
rough, order-of-magnitude estimates (for example, (Tr^r) 2a'^ ^^kr^)- For 
such states, then, we have an estimated ratio 



{\Sqi.Atf)\)^ 



Aqkritf) 

It then follows that if (instead of (HI])) 



Akaf^(ni,r)f + 1/2^ ^ dt ^V("kr) + 1/2 



1 / — ^ 

or equally if (instead of ([30 



j: ^ 4^af^{nur)f + 1/2^ ' dt ^^("kr) + 1/2 , (38) 



*/ . 1 



4afJaf{H^r)J dt -^^j a {Hkr) > I , (39) 



then 



(l^'Zk,.(i/)|) 



- > 1 • (40) 



Agkr(t/) 

From this we may reasonably deduce that relaxation, or at least significant 
relaxation, is likely to occur (except of course for special states with very simple 
velocity fields). 

Unlike our proof of relaxation suppression for modes satisfying this is 
not a rigorous result. (It is roughly analogous to saying, in classical kinetic 
theory, that significant relaxation to thermal equilibrium is likely to occur, over 
timescales of order the mean free time, if the mean magnitude of momentum 
transferred in molecular collisions is comparable to the width of the equilib- 
rium momentum distribution.) To delineate the precise behaviour in this region 
requires further study, perhaps through numerical simulations. 

To avoid potential misunderstandings, we should emphasise that relaxation 
might of course be suppressed for special quantum states violating the freezing 
inequality (|29|) (in particular, states with an especially simple de Broglie veloc- 
ity field). However, one should bear in mind that we are concerned with the 
evolution of quantum nonequilibrium in our actual universe, which is known to 
have had a complex and violent past history. Thus, for example, in a standard 
radiation-dominated phase, special states with no entanglement at any time are 
of no interest: we are concerned with states that are likely to have actually 
occurred. In seeking a general criterion for the freezing of early nonequilibrium, 
it is then of no use to point to special quantum states exhibiting particularly 
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simple velocity fieldsH In contrast, the freezing inequality is a natural 
constraint on quantum states in general, providing a realistic pointer to where 
nonequilibrium might be found in our actual universe. And violation of (j29p is, 
as we have argued in this section, likely to imply relaxation or at least significant 
relaxation. 



8 Time Evolution of (nk? ) for a Free Field 

For a given expansion history a = a{t) on [ti,tf], the freezing inequality (j29p 

depends on the time evolution of {nkr)t on [U, tf] (or, (PD|) depends on ^iJkr^ )■ 

To make precise predictions, then, we require a specific cosmological model, and 
an explicit expression for ( function of time t. We leave such detailed 

studies for future work. Here, we give a method for calculating (nkr)t (and 

Hkr^ ) for an arbitrary pure quantum state. This method might prove useful. 

The mean energy 



Wkr = (Hkr) = {{nkr) + l/2)(fc/a) 



in the mode kr evolves in time according to dW\i_r/dt — ^9-ffkr/9ty, which 
implies (using a = Ha) 

= -iHW^r + AHU^r , (41) 

where ?7kr = (s'^'^^^'^kr) mean potential energy. (For an interacting 

field, as discussed in section 5.2, dW^r/dt would contain additional terms from 

—i ([Hkr,Hi]^.) The rate of change of (rikr) = {a/k)Wkr — 1/2 is then given 

by 

= -2{Ha/k){K^r - f/k.) , (42) 

where i^kr = {'^''kr/'^^'^} mean kinetic energy. 

To solve for Wur{t) — K-kir{t) + C/kr(0; and hence the required function 
('"■kr)f , one may write first-order (linear) differential equations for JiTkr, f^kr and 

for the quantity Xkr = ^ ('Zkr'^kr + Ti'kr^kr)- Using d (Clj /dt = —i ([Cl, H]"^ 
dCl/dt^, it is readily shown that 

—— = -3HKkr-^ Xkr , —7— = HUkr+L^ Xkr , —J— = 2(Akr-(>kr-) • 

dt dt dt 

(43) 



•^A notable exception is the inflationary vacuum, which is in fact an example of a state 
that is non-entangled (across modes), with a very simple velocity field, and which is widely 
believed to have existed in the past; see section 10.1. 
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If H = a/a and uj — k/a are known functions of time, then given values of ii'kr, 
Uur, Xkr at any one time (say ti or tf) — where these values are determined 
by the wave functional 4" at that timqfl — the equations determine iiTkr, 
Ukr, Xkr at all times, yielding VFkr(i) = -fi'kr(i) + C^kr(i) as well as the required 
function (nkr)t = a{t)Wi^r{t)/k — 1/2. 

Introducing the vector X = (-ftTkr, C^kr, Xkr)"^, the equations take the 
form dX/dt = AX, where A is the time-dependent matrix 

/-3d/a -P/a^X 
A= a/a P/a^ . (44) 
\ 2 -2 / 

For interesting forms of a, such as a c>c i^/^, it seems likely that these equations 
will have to be solved numerically. 

It would be interesting to study this system of equations, and to establish 

the conditions under which solutions for (7^kr)f (or IVkr(^) = (ji'kr^ ) satisfy 

the freezing inequality (or ([50]) ). We leave this for future work. 



9 Approximate Solutions for (nkr)^ Satisfying the 
Freezing Inequality 

However, it is important to show first of all that solutions for {nkr)t satisfying 
(|29p can exist for some values of k. Here, we construct approximate solutions of 
(|43|l valid in the long- wavelength limit fc ^ 0, that satisfy ((29|l for appropriate 
initial conditions and time intervals. The conditions of validity are probably 
too restrictive for useful application to realistic cosmological scenarios, and we 
give these solutions here only to show that solutions satisfying are indeed 
possible. 

We consider a radiation-dominated expansion, for which a cx ^1/2 and H = 
l/2t. Dropping the indices kr, we find approximate solutions to (|43p satisfying 
(for appropriate values of k) 

w^lxl « HK, HU (45) 
(where K, U are non-negative), or 

^\X\«K,U (46) 

(where ti/a'^ — tf/a'j — tja^). We then have the simple solutions 

K^K.itlUY^I'^ , U = U,{t/uy/^ (47) 

"^Of course, initial values for -fCk^, f7kr, Xkr cannot be chosen completely arbitrarily. They 
are subject to constraints, such as a-^ K + U/{ak'^) + x = | (^{-k + g)^) > (or ak^K + U/{a^) + 

^^X > 0). 
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(that is, K (x. \/a^ and U (x a), and 

X = X^ + ^K,t, (l - + \u,U (l - {t/kf'^) ■ (48) 

Note that, for these solutions, the quantities (g^r) = '^U\^r/ (k^a) and (tt^^) = 
2a^iV'kr are time independent. 

We need to show the consistency of the solutions ((47l) and (|48)l with the 
assumed approximation (j46p . This may be done if k is appropriately small. 
Specifically, writing 

X^X^+ ^K,U + '^U^U - ^Kt - ^Ut , 
we have (since Kt and Ut respectively decrease and increase with time) 



IXI < \x^\+W^ + ^Uftf=D 



If we assume that 



k^ «-^^mm{Kf,U,} (49) 
(where (ti/D) min{A'/, Ui} is dimensionless), we then have 

^Ixl < «mm{Kf,U,} <K, U 

(since K and U respectively decrease and increase), and so the approximation 
condition ((46|) is indeed satisfied. 

For k satisfying (HS)) , we then have the approximate solutions (|T7l) for K and 
U. We wish to show explicitly that, for these solutions, there are values of k 
that satisfy the freezing inequality (I^S)) (or ([50)1 V 

To show this, for simplicity we first choose initial conditions with Ki << Ui. 
Since K decreases with time, we then have inm{Kf, Ui} — Kf and (from (|49l) ) 
the solutions ((T7|) are valid if 

« If (50) 

(where af/ti = aj^/tf — a? /t). Further, since K decreases and U increases with 
time, Ki « Ui imphes that K « U for all t > tj. Thus we have <^iJ^ « U{t) 
(where we continue to suppress the indices kr), or (using (|47|) ) 

h)^ ~ U,{t/uy/^ . (51) 

Inserting this into the freezing inequality ([30)) . and using a — af{t/tfy^^ and 
HJ^ = 2t-f, and taking ti/tf « 1, we obtain 

< 

4 af Hf 
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Since (H) Ui we have (restoring indices kr) the freezing inequaUty 



la, 1 1 / a, \ ^ 1 



(Note that, since for the above solution (^Hi^rj increases with time, the general 
result (I57|) also apphes, with (iJkr / = (Hkr) < This is consistent 

\ I i \ I min '' 

with (f53l) . since we have assumed tijtf « 1 which implies a,i/a/ << 1.) 

Thus, for a given mode kr satisfying ([5D1) . if is sufficiently small 

(satisfying (|53p ). then relaxation will be suppressed and initial nonequilibrium 
(if it exists) will be frozen. And it is indeed always possible to choose (^i?kr^ so 
as to satisfy ([55)1 . provided fc is sufficiently small. For the only general constraint 
on ^-ffkr^ is 

i?kr) = (("kr), + :r)- > ^- , 
/ I ^ Cii ^ Cii 



SO it is possible to satisfy ([53)1 if 



1 a, a. 
k < - 



2 a/ 
or 

Aphys = a/A > 47r i?^^ » HJ^ . 

If instead we choose initial conditions with Ki >> Ui, we will have K >> U 
only for as long as ti/tf is not much smaller than 1. Over this limited time, we 
have 

(h)^ « Kit) = miU)-^'^ = Kfit/tf)-^/^ . (54) 
Since Kf >> Uf > Ui, we now have min {-ft'/, Ui} = Ui and the solution is valid 

« f I . (55) 

Inserting ((54|) into the freezing inequality ([30]) . and assuming that U/tf is small 
compared to 1 (but not so small as to invalidate the approximation K >> U), 
we obtain the freezing inequality 

H^r)^^K..<\j^. (56) 

Again using (^Hkr'^ > (l/2)(fc/ai), we now find that it is possible to satisfy 
(ESI) if ' 

< 1 _ 1 ^a/ " 



2 i/ri 2 \a, J HJ^ 
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Aphys = a/A > 47r — H,. . 

af 

By assumption, ai/aj is not much smaller than 1, so we still have Aphys ^J^- 
(In any case, it follows from ([55)) that the solution is valid only if a/ A > HJ^. 
For we have 

8 8 
D > 8K^U + ^Uftf > 8K,U + ^U^U ~ SK,U , 

so that ([55]) gives 

U D tf 8K, \tj 
or afX >> {ai/af)HJ^ . Since aijaf is not much smaller than 1, we indeed have 

It is therefore certainly possible to have solutions for {nkr)t with nonequi- 
librium freezing for long-wavelength modes, a/A > HJ^ (or a/A >> HJ^). 



10 Possible Consequences of Early Nonequilib- 
rium Freezing 

The freezing inequality (|29p (or ([50)1 1 makes it possible, for the first time, to 
make quantitative predictions for nonequilibrium deviations from quantum the- 
ory, if we are given a specific cosmological model. The potential consequences 
are many, and much remains to be done to develop them. Here, we restrict our- 
selves to a preliminary sketch of some possible nonequilibrium effects, in par- 
ticular: corrections to inflationary predictions for the CMB, non-inflationary 
super-Hubble field correlations, and relic nonequilibrium particles. We hope 
to develop further details elsewhere, in the context of specific (and realistic) 
cosmological models. 

As we saw in section 6, for a radiation-dominated expansion (|29p implies the 
general lower bound ([55)) on the physical wavelength Aphys (i/) — the wavelength 
of what might be termed 'relic nonequilibrium field modes' — at the final time 
tf. In terms of the ambient temperature T, where T cx 1/a cx the lower 

bound may be written as 

Aphys(i/)>4^i/7'ln(I^./T/) . (57) 

As we have discussed, this lower bound will in practice be not more than two 
or three orders of magnitude larger than the Hubble radius H7^ at time tf. 

Note that, to satisfy the freezing inequality, the bound ([57[) is a necessary 
but not sufficient condition. A detailed understanding of where nonequilibrium 
freezing can occur requires, as discussed in section 5.1, a calculation of the time 
evolution of the mean occupation numbers (fikr)^ for the pure subensembles 
(with wave functionals ^'c) contained in the early mixed state, to find out which 
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— if any — of these subensembles satisfy ([^U]) . This is a matter for future work. 
Here, we consider only the necessary condition (j57p , which provides a pointer to 
where residual nonequilibrium could be found (pending the said more complete 
analysis). In particular, (j57p suggests that one should look for nonequilibrium 
above a specific critical wavelength. 



10.1 Corrections to Inflationary Predictions for the CMB 

In inflationary cosmology, the universe undergoes a period of exponential expan- 
sion, a{t) oc e^*, driven by the energy density of an approximately homogeneous 
scalar or inflaton field 0, where quantum fluctuations in (j) seed the primordial 
curvature perturbations that are later imprinted as temperature anisotropics in 
the CMB dH. 

To a first approximation, inflation predicts that modes of the inflaton field 
will have a quantum variance 

and a scale-invariant power spectrum 

where ( loikl^ / is obtained from the Bunch-Davies vacuum in de Sitter space, 
\ /qt 

for Aphys >> H^^. In the slow-roll limit {H 0), this results in a scale- 
invariant spectrum, {k) = const., for the primordial curvature perturbation 
7?.k, in approximate agreement with what is observed in the CMB [33) . 

Now, quantum nonequilibrium in the early Bunch-Davies vacuum generally 
implies deviations from ((58)l . It has been shown [lU [20] that if (microscopic) 
quantum nonequilibrium exists at the onset of inflation, then instead of relax- 
ing it will be preserved during the inflationary phase, and furthermore it will 
be transferred to macroscopic lengthscales by the expansion of physical wave- 
lengths Aphys oc a{t) cx e^*. Speciflcally, for each mode k, explicit calculation 
shows that the width of the evolving nonequilibrium distribution remains in a 
constant ratio with the width of the equilibrium distribution. (This is essentially 
because the vacuum state has the special property of being non-entangled across 
modes, so that the de Broglie-Bohm trajectories decompose into independent 
one-dimensional motions. See ref. [20].) If we write the nonequilibrium variance 
as 

(|0kP) = (|0kP)QT^(fc) (60) 

(where equilibrium corresponds of course to ^(fc) = 1 for all fc), the power spec- 
trum for 7?,k is then just the quantum result multiplied by the 'nonequilibrium 
factor' ^(fc): that is, Vn{k) = V^(k)i{k). 
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Thus, quantum nonequilibrium at the beginning of inflation will generally 
break the scale invariance of V-nik). As discussed in detail elsewhere [20], mea- 
surements of the angular power spectrum for the CMB may be used (in the 
context of inflation) to set bounds on ^(fc). 

Given these results, the next step is to try to predict some features of the 
function ^(fc). This requires a constraint on the form of nonequilibrium at the 
onset of inflation. 

One possible strategy is to consider a pre-inflationary era, and to derive con- 
straints on residual nonequilibrium from that era. If we take the pre-inflationary 
era to be radiation-dominated (a cx t^^'^), the lower bound ([57)) shows that 
nonequilibrium (for whatever fields may be present in that era) can survive 
only for sufficiently large, super-Hubble wavelengths. Since Aphys oc t^/^ and 
oc t, at sufficiently early times all physical wavelengths will in fact be 
super-Hubble (Aphys > H~^), raising the possibility of nonequilibrium freezing 
for the corresponding modes (if the freezing inequality ([29]) is satisfied). Dur- 
ing the subsequent inflationary phase, is (approximately) constant, and 
relevant cosmological fluctuations originate from inside . Some of these 
fluctuating modes could be out of equilibrium only if they evolved from modes 
that were outside the Hubble radius in the pre-inflationary phase. 

Thus, in order to obtain nonequilibrium corrections to inflationary predic- 
tions for the CMB, arising from an earlier pre-inflationary era, some of the pre- 
inflationary nonequilibrium modes must enter the Hubble radius, and they must 
avoid complete relaxation by the time inflation begins. Because pre-inflationary 
modes with larger values of A enter the Hubble radius later, they are presumably 
less likely to relax before inflation begins, in which case residual nonequilibrium 
will be possible only for A larger than some infra-red cutoff Ac. (For further 
discussion, see ref. [20] ■) 

We hope that future work, based on a specific pre-infiationary model, will 
provide a prediction for Ac, as well as some indication of the form of the nonequi- 
librium spectrum for A > Ac. Note that ^(/c) < 1 at wave number k implies that 
the nonequilibrium width of the corresponding inflaton mode is less than the 
equilibrium width. One might reasonably expect this, in view of the hypothesis 
that quantum noise arose from statistical relaxation processes in the very early 
universe: it seems natural to assume that early nonequilibrium would have a 
less-than-quantum dispersion, ^(fc) < 1, as opposed to a larger-than-quantum 
dispersion, ^(k) > 1 (though the latter is of course possible in principle). Thus, a 
dip f (fc) < 1 in the power spectrum below some critical wave number kc = 27r/Ac 
might be naturally explained in terms of quantum nonequilibrium surviving from 
a very early pre-inflationary era. 

It has in fact been found that an infra-red cutoff in the primordial power 
spectrum provides a slightly better fit to the 3-year WMAP data; however, 
the improvement is not sufficient to justify introducing the additional cutoff 
parameter in the model [34] . 
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10.2 Super-Hubble Correlations without Inflation? 

As noted in the introduction, one motivation for assuming quantum nonequi- 
librium at the big bang was that the resulting nonlocahty at early times could 
eliminate the cosmological horizon problem (which persists, as we have men- 
tioned, even in some inflationary models [IH])- One might also ask if early 
quantum noncquilibrium could provide an alternative, non-inflationary means 
of laying down primordial curvature perturbations at super-Hubble lengthscales, 
in a standard Friedmann cosmology. Since we have shown that noncquilibrium 
can remain frozen at super-Hubble scales, one may ask if such noncquilibrium 
could generate appropriate super-Hubble correlations without the need for an 
inflationary era. 

The Bunch-Davies vacuum for a scalar field (j>, with variance given by ([58)) (at 
long wavelengths), has the remarkable property that the two-point correlation 
function 

(0(xi)0(x2)) = ly d^k e'-^-^-^--^' (|0kP) 

is independent of distance |xi — X2I, as is readily verified for (l^kP) = (l0kP)Qrp oc 

As a first step, one may ask how this inflationary quantum behaviour 
could be mimicked by a non-inflationary vacuum in quantum noncquilibrium. 

Consider a vacuum state whose quantum variance is (|</'kP)Qrp oc fc™^'^ 
for some fixed index toqt. Assuming that the quantum two-point function 
decreases with distance, where ((/)(xi)0(x2))Qrp oc |xi — X2I ^'"'3'^+'^)^ we have 
?7iQT > —3. Then consider the same vacuum state in quantum nonequilibrium, 
with (|(?I)kP) = (l0kP)Qrp assuming that ^{k) oc for some fixed index /i. 
To obtain a nonequilibrium two-point function that is independent of distance, 
we require toqt + ^jl — —3, or < 0, so that (in this simple example) the 
nonequilibrium function ^(fc) must increase as fc ^ 0. 

As things stand, we are unable to say if such behaviour for £^{k) is likely to 
emerge from any reasonable model. However, given the upper bound (j28p on the 
ratio (|(5qkr(i/)|)cq /Aqkr(i/), one could study how the 'degree of freezing' varies 
with k (for example for A; ^ 0), where a high or low degree of freezing could 
be defined respectively as a low or high value of the upper bound on the right- 
hand-side of For a specific cosmological model, with some assumptions 
about initial conditions, this could provide constraints on the behaviour of the 
function ^(fc). The results will obviously depend on how (rikr) varies with k. 

Finally, we note that a nonlocal model, based not on hidden variables or 
quantum nonequilibrium but on the holographic principle, has been shown to 
generate the required (approximately) scale-invariant perturbation spectrum at 
super-Hubble scales |35| . Whether or not early quantum nonequilibrium could 
reproduce such effects in a natural way remains to be seen. 
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10.3 Relic Nonequilibrium Particles 

We saw in section 10.1 that relic nonequilibrium field modes could, in the case 
of the inflaton, change the power spectrum for primordial curvature perturba- 
tions, resulting in observable effects in the CMB. Thus, inflationary cosmology 
provides a simple and definite means whereby early nonequilibrium could yield 
observable consequences today. In this section, in contrast, we shall attempt 
to outline some much more complicated and uncertain scenarios, according to 
which relic nonequilibrium field modes (for some appropriate field) might in 
some circumstances manifest as relic nonequilibrium particles that could be de- 
tected today. Unfortunately, for these scenarios to be at all plausible, some 
questionable assumptions have to be made, and at the time of writing it is not 
clear if these scenarios can really work in practice. 

There is of course no preferred definition of particle states in quantum field 
theory on expanding space, except in the short-wavelength limit (where one 
recovers the usual Minkowski definition) [35]. The very notion of 'particles' is 
in fact highly ambiguous for modes of frequency lower than the typical inverse 
timescale over which the spacetime metric changes. In a cosmological setting, 
this means that there is no generally useful definition of quantum particle states 
at wavelengths larger than the Hubble radius. Thus, if we consider relic nonequi- 
librium field modes from a radiation-dominated era — where the bound (|57p 
implies that such modes must have super-Hubble wavelengths — we must be 
careful not to interpret such modes too naively in terms of (quantum) particle 
states. However, pending a more precise treatment, one might reasonably as- 
sume that if such modes enter the Hubble radius at later times, then they will 
manifest as (approximately-defined) particle states in the usual sense. 

One should also bear in mind that, generally speaking, excitations of super- 
Hubble modes will not be produced by the local processes of particle scattering 
and decay (which are not expected to be effective over lengthscales larger than 
the instantaneous Hubble radius H~^). However, such excitations will of course 
be produced by the global effects of spatial expansion. 

In order to maximise the chance of obtaining relic nonequilibrium particles 
that could be detected in practice (in particular, with energies that are not so 
low as to be completely out of range), we ought to try to minimise the lower 
bound on the mode wavelength defined by ([35]) or ([57]) . This can be done by 
choosing the final time tf to be as small as possible, subject to the constraint 
that further relaxation may be neglected for times later than the chosen value 
oitf. Thus, one might take tf to be the time tdcc at which the relevant particle 
species decouples. For one might reasonably assume that relaxation may be 
neglected (at all wavelengths) for t > tdec — if the quantum states, defined post- 
decoherence, are such that the associated de Broglie velocity field is sufficiently 
simple (as occurs, for example, for energy eigenstates). For a super-Hubble 
mode at idee that becomes sub-Hubble at later times, it is then conceivable 
that any nonequilibrium present at the time tdcc could persist until much later. 
(A proper discussion of this scenario would require an analysis of decoherence 
before and after decoupling.) 
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If we make the above assumptions, the key question is then whether resid- 
ual nonequihbrium field modes can exist at the time tj ~ idoc- From (|57[) . 
this is possible if the modes have physical wavelength (inserting the Boltzmann 
constant /cb) 

Aphys(idoc) > 47ri/^^^ln(fcBTi/A:BT'dcc) , (61) 

where H~[^^ and T^ec are respectively the Hubble radius and temperature at 
time tdcc- We have Aphys(idec) = adccA, where Odcc = To/Tdec (with Tq ~ 2.7 K 
the temperature today). Assuming that decoupling occurs before the end of 
the radiation-dominated phase, we also have H^^^ = 2tdec, where fdcc may be 
expressed in terms of Tdec using the standard temperature-time relation 



ksT J 

The lower bound ([6T|) then becomes (inserting the speed of light c) 

A>8.c(ls)fi^Vi^VnfT^V (63) 

\kBTdccJ \ Kb Jo / VKBTdec/ 

or (withe ~ 3x10^° cms^i, /cb ~ 8.6x10"'^ eVK^\ and /cbTq ~ 2.3x10"^ eV) 

This provides a lower bound on the wavelength A today, at which noncquilibrium 
could be found. 

The freezing inequality ((29|) . and the resulting lower bound ([64]), have been 
derived in this paper for massless scalar fields only. One certainly expects to 
find comparable results for more general massless boson fields, such as the elec- 
tromagnetic field. For fermions, however, a separate analysis is required. There 
are different approaches to the pilot-wave theory of fermions, and the details 
of nonequihbrium freezing may depend on which model is adopted. One might 
try to derive a fermionic analogue of the freezing inequality using, for example, 
the Dirac sea pilot-wave model [37|. Pending such extensions of our analysis, 
here we assume that the lower bound (|64p applies (at least approximately) to 
fermions as well, provided they are effectively massless at the temperature Tdcc 
(that is, of mass m « fcB?doc/c^)- 

With this understanding, let us now apply the approximate result (|64[) to 
various particle species, both bosonic and fermionic. For definiteness, we first 
consider a standard Friedmann cosmology with no inflationary period, taking 
our initial conditions at the Planck era /cbT^ ^ A:bTp ^ 10^^ GeV. (An alter- 
native possibility, of nonequihbrium relic particles arising from the decay of the 
infiaton, is considered below.) 

Photons decouple from matter at fcB(Tdcc)7 ^ 0.3 eV. From ([64]) we then 
have a lower bound 

A.^ > 0.7 X 10^° cm , (65) 
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which exceeds the Hubble radius today, Hq ~ 10 cm. If instead we consider 
neutrinos, which decouple at fcB(T'doc)i/ ^ 1 MeV, we have 

> 1.7 X 10^3 cm ~ 5.5 x lO"^ pc (66) 

(or ~ 10^ light years). Residual nonequilibrium for relic neutrinos could plau- 
sibly exist today only at such tiny energies. Unfortunately, this is of course far 
outside any realistic range of detection. (Note, again, the implicit assumption 
being made, that if nonequilibrium super-Hubble modes at tdec enter the Hubble 
radius at < > idee, they will manifest as nonequilibrium particle states.) 

The situation improves drastically, however, if one considers particles that 
decouple soon after the Planck era. Gravitons, for example, are expected to 
decouple at a temperature (rdcc)g ^ 7p. Writing 

fcB(rdcc)<, = XgiknTp) ~ XgilO'^ GeV) , 

where Xg < 1, we obtain 

Ag> (0.3cm)(l/x<,)ln(l/a;g) . (67) 

This might be compared with the range of wavelengths expected for a (ther- 
mal) relic graviton background, whose temperature today is estimated to be 
(2o)g ~ 1 K [35j. At this temperature, the spectral energy density of a Planck 
distribution peaks at the wavelength A„iax(l K) ~ 0.3 cm. 

There may also exist other particles that decouple not too long after the 
Planck era, and that (unlike the graviton) are unstable, eventually producing 
decay products that could be more easily detected today. A natural candidate, 
arising out of current supersymmetric theories of high-energy physics, is the 
unstable gravitino G, which has been estimated to decouple at a temperature 




where is the number of spin degrees of freedom (for the effectively massless 
particles) at the temperature (rdcc)(=;, nT-gi is the gluino mass, and to^ is the 
gravitino mass. This provides us with an estimate for the lower bound in the 
case of gravitinos, 

AG>(0.3cm)(l/x^)ln(l/a;^) . (68) 

For the purposes of illustration, if we take ((7*/230)^^^ ^ 1 and (1 TeV /mgif' ~ 
1, then 

~ Vl03 GeV J ■ 

If, for example, tuq « 100 GeV, then Xq ~ 10~^ and ([SS]) yields Xq > 140 cm. 
This corresponds to energies that are rather low, but perhaps accessible. 

If the gravitino is not the lightest supersymmetric particle, then it will indeed 
be unstable. For large m^, the total decay rate is estimated to be [40] Tq = 
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(193/48)(m|/M^), where Mp ~ 1.2 x 10^^ GeV is the Planck mass. The time 
(idocay)(5 at which the gravitino decays is of order the hfctime l/r^. Using ((62l) . 
the corresponding temperature is 

fcB(rdecay)G ~ GeV)^/^ eV . 

For example, again for the case m^, « 100 GeV, the relic gravitinos decay 
when fcB(Tdecay)(5 ^ 1 kcV. This is prior to photon decoupling, so that any 
(potentially nonequilibrium) photons produced by the decaying gravitinos would 
interact strongly with matter and quickly relax to quantum equilibrium. To 
obtain gravitino decay after photon decoupling, we would need kB{Tdccay)Q ^ 
kBiTdeclj ~ 0.3 eV, or m^, < 0.5 GeV. For such small gravitino masses, 
however, decoupling occurs at (roughly) 

(Tdoc)G = x^Tp w (tog/IO^ GeV)'Tp < IQ-^Tp 

and (with Xq < 10"'') yields the much larger lower bound Xq > 10'' cm. 
Thus, it may prove more promising to consider other decay products (that 
decouple prior to gravitino decay but for larger gravitino masses). These could 
in turn decay into photons at later times, or they might be detected directly. 

There are of course strong constraints on the presence of gravitinos in cos- 
mological models, in particular from the abundance of light elements emerging 
from big-bang nucleosynthesis and from limits on dark matter abundance. These 
constraints have been extensively studied — see, for example, ref. [41j — and 
the subject is an active area of current research. Our hope is that an accept- 
able and compelling scenario will eventually be found, satisfying the standard 
cosmological constraints and at the same time allowing the possibility of relic 
nonequilibrium surviving in particles that could be detected today. To develop 
such a scenario in detail is a topic for future work. 

So far in this section, we have assumed a standard (non-inflationary) Fried- 
mann expansion, with initial nonequilibrium at around the Planck era. An 
alternative scenario is obtained if we consider relic nonequilibrium particles in 
the context of inflationary cosmology. If inflation did occur, the density of any 
relic particles (nonequilibrium or otherwise) from a pre-inflationary era will of 
course be so diluted as to be completely undetectable today. However, one may 
consider relic particles that were created at the end of inflation, by the decay of 
the inflaton field itself. 

As discussed in section 10.1, during inflation the inflaton field does not relax 
to quantum equilibrium, and in fact the exponential expansion of space transfers 
any initial nonequilibrium from microscopic to macroscopic lengthscales. The 
inflaton field, then, is a prime candidate for a carrier of primordial quantum 
nonequilibrium. As well as manifesting as statistical anomalies in the CMB, 
such nonequilibrium in the inflaton field could manifest as nonequilibrium in 
its decay products, where in standard inflationary scenarios inflaton decay is in 
fact the source of the matter and radiation present in our universe today. 

The process of 'preheating' is driven by the homogeneous and essentially 
classical part of the inflaton field (that is, by the fc = mode) [42]. Here, the 
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inflaton is treated as a classical external field, acting on other (quantum) fields 
which become excited by parametric resonance. Because of the classicality of 
the relevant part of the inflaton field, this process is unlikely to result in a 
transference of nonequilibrium from the inflaton to the created particles. 

During 'reheating', however, perturbative decay of the inflaton can occur, 
and one may reasonably expect nonequilibrium in the inflaton field to be trans- 
ferred to its decay products. This possibility opens up a large field of investiga- 
tion. Here, again, we restrict ourselves to making some preliminary remarks. 

The perturbative decay of the inflaton occurs through local field-theoretical 
interactions, so one expects the decay products to have physical wavelengths no 
greater than the instantaneous Hubble radius. Taking the lower bound ((57|) as 
a guide (even though it was derived for a radiation-dominated phase), we then 
expect that the decay products will come into existence already violating the 
freezing inequality. Subsequent relaxation might then be avoided (possibly) only 
if the particles are created at a temperature below their decoupling temperature. 
Once again, the gravitino suggests itself as a possible candidate. Gravitinos can 
in fact be copiously produced by inflaton decay (and could even make up a 
significant component of dark matter [44]). If the gravitinos are unstable, again, 
one could try to detect (say) photons produced by their decay at later times. 

The possible realisation of this scenario depends of course on uncertain fea- 
tures of high-energy particle physics and of inflationary models. As before, one 
may hope that a scenario will eventually be found, satisfying the constraints of 
particle physics and cosmology, and at the same time allowing the possibility of 
relic nonequilibrium surviving in particles that could be detected today. 

We close this section with some general remarks. 

First, we note that particle decay (for example for the gravitino) is likely to 
result in some relaxation and erasure of any quantum nonequilibrium that may 
have existed in the parent particles. However, one hopes that the erasure will 
not be complete and that some nonequilibrium will still be present in the decay 
products. It would be useful to study this, in pilot- wave models of specific decay 
processes. 

Second, once suitable candidates for nonequilibrium relic particles have been 
identified, one must consider how best to test them for violations of the Born 
rule. For photons, a particularly simple test involves searching for anomalous 
polarisation probabilities, or deviations from Malus' law (where such devia- 
tions reflect the nonequilibrium breakdown of expectation additivity for non- 
commuting quantum observables in a two-state system) [121 113 • 

Third, for a given species of relic particle in the universe today, even if there 
exist pure subensembles with significant residual nonequilibrium, in practice it 
might be difficult for us to locate those subensembles and perform experiments 
with them. In particular, if a given detector registers particles belonging to 
different subensembles, without distinguishing between them, it is possible that 
even if nonequilibrium is present in the individual subensembles it will not be 
visible in the data. 
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11 Conclusion 



The hypothesis of quantum nonequiUbrium at the big bang has been shown 
to have a number of observable consequences. Our main resuh is the freezing 
inequahty For cosmological field modes satisfying initial nonequi- 

librium will be 'frozen' at later times. This result may be applied to specific 
cosmological models, yielding predictions whose verification could constitute ev- 
idence for quantum nonequilibrium in our universe. For a radiation-dominated 
expansion, ((29|) implies the general lower bound p5|) on the wavelength of relic 
nonequilibrium field modes. 

The detailed study of quantum nonequilibrium freezing, for realistic cosmo- 
logical models, is left for future work. A useful first step might be to study the 
system of equations (j43p . and to delineate the general conditions under which 
the time evolution of a (mean) mode occupation number {nkr)t can satisfy the 
freezing inequality (P^ . Crucially, future work will need to study the statistical 
distribution of wave functionals for a realistic mixed state on expanding space, 
the goal being to identify subensembles satisfying (P^ . For these subensembles, 
quantum nonequilibrium is expected to be frozen over the relevant time period, 
resulting in definite predictions that might be tested today. 
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